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ABSTRACT 

The  specific  heat  critical  exponent  a  has  been  measured  for  the  A u(  1 10)  (1x2) 
order-disorder  phase  transition  using  partially  integrated  Low-Energy  Electron 
Diffraction  (LEED)  intensities.  The  resulting  value,  a  -  0.02  ±  0.05,  Is  consistent 
with  the  predicted  Islng  universality  class  of  this  transition.  Evidence  for  a 
reduction  in  the  effective  critical  temperature  due  to  finite  size  effects  Is  also 
presented. 
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Studies  of  phase  transitions  on  real  surfaces  provide  Important  tests  of  our 
basic  Ideas  about  two  dimensional  critical  phenomena.  The  Au(  1 1  OK  1  x2)  surface  Is 
an  excellent  model  for  studies  of  continuous  order-disorder  transitions.  The  clean 
surface  Is  inert  to  the  residual  gases  normally  present  in  an  ultra-high  vacuum 
chamber  so  that  much  longer  times  are  available  for  measurements  than  is  usually 
the  case.  The  phase  transition  is  known  to  belong  to  the  Ising  universality  class  and 
the  critical  exponents,  13,  y  and  v  have  been  measured  by  Campuzano  et  al.  [  I  ].  Thus, 
Au(  1 1  OK  1  x2)  can  be  used  for  more  ref ined  studies  of  critical  behavior  including  the 
role  of  finite  size  effects  [2]  and  the  kinetics  of  ordering  [3],  In  this  report  we  use 
Au(  1 1 0)(  1  x2)  to  verify  a  method  proposed  by  Bartelt,  Einstein  and  Roelofs  14]  for 
extracting  the  heat  capacity  critical  exponent  a  from  low-energy  electron  diffraction 
(LEED)  data.  We  also  find  evidence  for  finite  size  effects  In  this  system.  It  is  also 
significant  that  the  method  of  BER  significantly  expands  the  ability  of  LEED  to  study 
phase  transitions  since  the  required  measurements  and  analysis  are  simple  and 
multiple  scattering  is  fully  accounted  for. 

BER  consider  LEED  intensities  Integrated  over  a  small  region  near  the  center  of 
the  surface  Brillouln  zone.  This  integration  reduces  the  sensitivity  to  long  range 
order  so  that  the  singular  behavior  near  a  phase  transition  Is  that  of  the  average 
energy  E.  Thus,  the  exponent  a,  defined  near  the  critical  temperature  Tc  by  [4] 

1-cu 

E-E0*a,t±bj|t|  '  ♦•••  (1) 

where  E0,  a|,  b4  and  b_  are  constant  coefficients,  t  ■  (T-Tc)/Tc  is  the  reduced 
temperature  and  ♦  and  -  refer  to  t  >  0  and  t  <  0  respectively.  Scaling  predicts  that  a ♦ 
■  a_  and  the  ratio  b4/b_  Is  a  universal  quantity.  For  the  Ising  universality  class,  a*  * 
0  which  implies  b4/b_  -  1.  As  pointed  out  by  BER,  most  real  LEED  instruments  have  a 
finite  transfer  width  to  and  are  only  sensitive  to  correlations  over  a  finite  range  Lj  “ 
co  so  that  the  required  average  is  obtained  automatically.  Furthermore,  since 
multiple  scattering  Is  short  ranged  [4],  Its  presence  in  LEED  does  not  hinder  the 


analysis.  BER  show  explicitly,  that  as  long  as  the  diffracted  intensity  contains  no 
Information  about  the  phase  of  the  order  parameter,  the  measured  LEED  Intensity  has 
the  form 

,  %  Aq-  Ajt  ♦  B.ltp-a-  ♦  t  <  0 

l(t)  =  (2) 

Aq  -  Ajt  -  B+ |t|'~a+  ♦  -  t>0. 

where  Aq,  Aj,  B+  and  B.  are  constants.  As  for  Equation  ( 1 ),  a4  ■  and  B+/B.  = 
b+/b_. 

BER  used  Monte  Carlo  simulations  for  (</3xV3)R30  and  p(2x2)  ordered  phases  on  a 
triangular  lattice  to  show  that,  depending  on  the  range  of  t.  Integration  of  the 
intensity  over  more  than  about  0.83%  of  the  Brillouin  zone  were  sufficient  for 
Equation  (2)  to  be  valid.  Since  the  integral  of  the  intensity  over  the  entire  Brillouin 
zone  is  conserved  [5],  the  radius  of  the  integration  must  also  be  much  smaller  than 
the  diameter  of  the  zone. 

The  largest  value  of  t  for  which  Equation  (2)  is  valid  Is  determined  by  the 
instrument  and  the  influence  of  corrections  to  scaling  terms;  the  smallest  depends  on 
the  perfection  of  the  surface  lattice.  The  smallest  Itl  =  tmjn  is  obtained  when  the 
correlation  length  t  becomes  comparable  to  the  length  scale  L  of  the  finite  size 
regions; 

(3) 

where  is  the  order  of  several  lattice  spacings.  A  more  perfect  surface  results  In  a 
small  tmjn  The  largest  value  Itl  *  tmax  occurs  when  i  is  comparable  to  the  transfer 
width  of  the  instrument,  ©; 

tmax“<“/<o>',/v.  (4) 

unless  there  are  important  corrections  to  scaling  beyond  those  indicated  in  Eq.  (I). 


The  Au(  1 10)  sample  is  the  same  one  used  in  a  Rutherford  Backscattering  study 
[6],  It  was  oriented  to  within  ±0.5’  of  the  (110)  plane  and  cleaned  in  ultra-high 
vacuum  with  1  keV  argon  ion  sputtering  followed  by  annealing  at  1070  K. 

Temperature  was  measured  with  two  type  K  thermocouples  [7]  held  against  the  side 
of  the  crystal  by  Pt  mounting  wires.  The  temperature  was  maintained  constant  to 
within  ±0.2  K  during  each  measurement.  Diffracted  intensities  were  measured  with  a 
movable  Faraday  collector  whose  circular  aperture  subtended  3.0x1 0‘5  sterradians  at 
the  sample.  A  retarding  grid  behind  the  aperture  was  biased  0.9V  below  the  primary 
beam  voltage. 

Figure  One  shows  the  intensity  of  the  (0,-fc  )  beam  as  a  function  of  temperature 
for  a  62  eV  primary  beam  incident  at  45’  along  the  [1 10]  direction.  The  inflection 
point  provides  an  estimate  of  Tc  «  695K  [4],  Below  650K,  the  measured  full-width  at 
half-maximum  (FWHT1)  of  this  beam  was  0.43  ±  0.02  nrn-1  in  good  agreement  with  the 
value  0.42  nm"1  estimated  for  our  LEED  instrument  using  the  procedure  described  by 
Park  et  al.  [8],  For  the  results  presented  in  Figures  1  and  2,  this  Instrument  response 
corresponds  to  integration  over  a  range  of  scattered  wave  vectors  equal  to  about  5.6% 
of  the  (1x2)  Brillouin  zone  length  along  the  [110]  direction.  Ratios  of  the  FY/HI 
measured  for  the  (0,-^  )  and  (0,- 1 )  beams  were  also  In  good  agreement  with 
calculated  values  and  provide  evidence  that  the  beam  width  was  dominated  by 
instrumental  effects  rather  than  finite  size  effects.  This  was  not  the  case  for  Ref.  I. 

Figure  One  also  shows  examples  of  a  (0,-!£ )  beam  angular  profile  measured  at 
478  K  by  fixing  the  angle  of  Incidence  and  the  detector  position  and  varying  the 
incidence  energy.  Possible  artifacts  associated  with  this  method  are  described 
elsewhere  /9/  but  are  not  Important  for  the  analysis  reported  here.  The  intensity  far 
away  from  the  maximum  was  used  to  define  a  linear  background  Intensity  due  to 
thermal  diffuse  scattering  (TD5)  and  Intrinsic  surface  imperfections.  Profiles 
measured  at  temperatures  far  enough  above  Tc  so  that  fluctuations  in  the  (1x2)  phase 
no  longer  contribute  significantly,  show  this  linear  approximation  to  be  accurate  to 
better  than  10*  of  the  background  Intensity  After  background  subtraction,  the 
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exponential  Debye-Waller  dependence  of  the  intensity  was  determined  using  data  with 
T  <  600  K.  The  effective  Debye  temperature  was  found  to  be  102.5  K  in  good 
agreement  with  previous  estimates  [6,10],  This  Debye  temperature  dependence  was 
assumed  to  apply  over  the  entire  temperature  range  of  the  data,  T  <900  K. 

We  used  the  following  procedure  to  analyze  the  critical  behavior  of  the  corrected 
data.  Tc  was  systematically  varied  between  640  K  and  740  K.  For  each  choice  of  Tc, 
Aq  was  chosen  by  a  linear  extrapolation  between  the  nearest  data  points  on  each  side 
of  Tc.  The  data  were  then  plotted  on  a  In  (KO-Aq)  vs  In  Itl  graph  and  the  longest 
straight  line  portion  on  each  side  of  Tc  was  determined  by  a  linear  least-squares 
fitting.  Examples  for  Tc  values  of  680  K,  695  K,  and  710  K  are  shown  in  Figure  Two. 
The  best  overall  fit  was  taken  to  be  that  which  most  closely  yielded  a+  =  a_  and  had 
nearly  equal  correlation  coefficients  for  t  >  0  and  t  <  0. 

The  best  fits  were  obtained  for  692  K  <  Tc  <  698  K  and  the  average  values  of  cu 
and  a_  over  this  range  of  Tc  was  a  *  0.02  t  0.05.  This  result  Is  in  excellent 
agreement  with  the  expected  Islng  value,  a  »  0. 

We  also  estimated  the  critical  exponent  v  =  1.1  ±  0. 1  from  the  measured  increase 
In  the  FWHM  of  the  (0,-*£ )  beam  In  the  range  0.02  <  t  <  0.07.  This  result  agrees  with 
the  Islng  value,  v  =  I  and  gives  added  confidence  in  our  results  for  a  and  Tc. 

Good  fits  for  a  were  obtained  for 
0.004  <  t  <  0.035. 

Mr  Mr 

The  diameter  of  finite  size  regions,  as  estimated  from  Equation  (3),  is  therefore 
about  250  Au(  1 1 0)  lattice  spacings  along  1 1 1 0]  or  about  1 00  nm.  This  result 
indicates  that  the  substrate  was  of  very  good  quality  and  that  its  contribution  to  the 
measured  beam  FWHM  is  negligible  compared  to  the  instrumental  broadening.  The 
value  of  a  is  not  very  sensitive  to  tmax  as  also  pointed  out  by  BER. 

We  note  that  the  analysis  described  above  has  in  effect  neglected  the  linear  term 
Ajt  In  Equation  (2).  This  is  necessitated  by  the  fact  that  a  is  zero  in  the  Islng  case 
so  that  the  critical  singularity  is  also  linear  In  Itl  to  lowest  order.  The  results  of  BER 
suggest  that  this  approximation  causes  the  value  of  a  extracted  from  the  analysis  to 
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be  too  large.  Although  our  analysis  does  not  enable  the  coefficients  B*.  to  be 
extracted  Inspection  of  Figure  2  shows  that  (A|  ♦  B+)/Aj  ♦  B_)  Is  nearly  unity.  A 
more  sophisticated  fitting  procedure  Is  not  warranted  given  the  accuracy  of  the  data. 

The  critical  temperature  of  695  t  3  K  is  substantially  larger  than  the  650  t  1.5  K 
obtained  by  Campuzano  et  al.  [1],  This  difference  is  consistent  with  that  expected 
from  finite  size  effects.  The  effective  (or  pseudo)  critical  temperature  Te  Is 
expected  to  vary  with  the  scale  L  of  finite  size  regions  approximately  as  [1 1] 

Tc  -  Te  ■  AT  «  aL"^  Tc. 

where  a  Is  a  constant  and  L  is  measured  in  lattice  spacings  (0.408  nm  along  the  [110] 
direction  for  Au(  1 10)). 

Landau  [12]  has  determined  a  for  the  case  of  a  finite  simple  Ising  model  with 
free  boundary  conditions  using  Monte  Carlo  simulations  and  finds  a  =  1.25  ±  0.04 
verifying  an  earlier  result  of  Ferdinand  and  Fisher  [1 1]  who  also  present  general 
arguments  that  a  should  be  of  order  unity  independent  of  the  details  of  the 
Hamiltonians  Assuming  that  our  measured  Tc  is  close  to  that  for  an  ideal  surface  and 
taking  Te  =  650  K  from  Campuzano  et  al.  { 1  ]  yields  L  “  20  lattice  spacings.  This  is  In 
excellent  agreement  with  the  mean-size  estimated  for  ordered  (1x2)  regions  on  the 
Au  crystal  used  by  Campuzano  et  al.  [1],  They  observed  the  full-width  at  half¬ 
maximum  of  the  fractional  order  diffraction  beam  to  be  about  0.05  of  the  spacing 
between  integral  order  beams  which  also  corresponds  to  ordered  regions  about  20 
substrate  lattice  spacings  wide. 

In  conclusion,  we  have  measured  the  specific  heat  exponent  for  Au(l  10X1x2) 
using  the  method  of  Integrated  intensities  proposed  by  Bartelt,  Einstein  and  Roelofs 
[4],  The  resulting  a  *  0.02  ♦  0.05  is  in  excellent  agreement  with  the  predicted  values 
and  thus  verifies  the  method  of  Integrated  Intensities  This  result  also  Illustrates 
that  the  method  greatly  Increases  the  utility  of  LEED  as  a  technique  to  study  critical 
phenomena  we  have  also  presented  evidence  for  a  dependence  of  the  effective 
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critical  temperature  due  to  finite  size  effects  consistent  with  theoretical 
predictions. 
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Figure  One:  (a)  Measured  Intensity  vs.  Temperature  before  (o)  and  after  (+) 
correction  for  the  Debye-Waller  factor,  (b)  (0,-^  )  Diffraction  beam  profile  at  T  - 
478  K.  The  background  correction  used  in  the  analysis  is  indicated  by  the  dashed  line. 

Figure  Two:  Plots  of  corrected  intensities  versus  reduced  temperature.  The  dashed 
lines  show  the  behavior  expected  for  an  Ising  system  with  a  =  0  and  the  solid  lines 
the  behavior  for  a  three  state  Potts  system  with  a  *  0.33. 


